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ABSTRACT 


This initial study investigates the feasibility of using 
weather map types to forecast tropical storm/typhoon movement 
in the Western North Pacific Ocean. Fleet Numerical Oceanog- 
raphic Center's hemispheric D-value analysis fields at 850, 
700 and 500 mb are interpolated to a standard grid oriented 
relative to the present surface tropical storm position. Enm- 
Dirical Orthogonal Function (EOF) analysis is used to repre- 
See @esrlelas in terms of only 10 coefficients in an effort 
to separate true synoptic map features from random noise. A 
least squares approach 1s then applied to determine character- 
imeele SyYnOptic patterns relative to the tropical storm, and 
to relate each individual case to a "map type". An analog 
type approach is used to forecast the 30-, 36-, 54-, 650-, 78- 
and 84-hour positions. All storms determined to be of the 
Same map type as the candidate storm were considered to be 
analogs, and their storm tracks were rotated to have the same 
past i2-nour movement. Forecast errors ranged from 135 n.mi. 
emeeeoours tO 490 n.mi. at 84 hours. Although no indepen- 
dent cases were tested, and the techniques employed have not 
been optimized, the mean vector errors indicate that the basis 


technique warrants further investigation. 
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The standard grid used to plot the D-value 
f1iéids. The standard position of the storm 
on every map 18S at position 70 - ---+e-+-+-+-+ - 12 


The forecasts 1, 2, 3 and 4 represent the 4 
scheduled warnings issued by JTWC each day. 

The X marks upper level analysis time; o 

represents time that forecasts would be 

available from this research, * represents 

the time warnings are issued; and $ repre- 

sents warning forecast times. This diagram 
1llustrates the reason for the 30-, 36-, 54-, 

60-, 78- and 84-hour forecast intervals - - - - - 1l/ 


This 1s a very generalized picture of the 

flow as it appears on the grid. Note the 

difference in speed, relative to the storm, 

of features to the north of the storm com- 

pared to features south of the storm - - --- - 18 


The effect of the different normalization 

schemes at 500 mb. The four representa- 

meomceware cor O00 GMT, 18 Mar 1967. A is 

mieerdernd! field; B, C and D are the EOF 
representations (truncated at 10 terms) 

using grid-point normalization, whole field 
normalization and map normalization, respec- 

tively. The isopleths are D-values (m) - - - - - 20 


Similar to Fig. 4 except for 12 GMT 
23 how 1 ee 2 ee ee 


The D-value means (left side) and standard 
@eweiataoens (right side) for 850 mb (bottom) , 

Poteme (middie) and 500 mb (top). The 

isopleths are D-values (m) - -----+-+-+--- 24 


An example of the effect of truncation of 

imicmees ror O00 GMT, 18 Mar 1967. The 

number above the map is the number of terms 

imelagded (ithe truncation point). The 

isopleths are D-values (m) - ------+--- - 30 
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Neumann (1979) states that the temporal and spatial 
enalysis of tropical cyclone paths show a tendency to be re- 
petitive, and to be associated with identifiable weather 
aee-ris. Experienced tropical cyclone forecasters routine- 
meexamine the synoptic fields for the position of the sub- 
maepteal ridge, any weakness in the intensity of the ridge, 
locations of migratory troughs in the mid-latitudes and any 
other synoptic features that may give indications of future 
storm movement. Ina series of lectures at the Naval Post- 
Beaduate School, Dr. William M. Gray, Colorado State Univer- 
Sity, indicated that much more can be done to examine how 
tropical storms react to their environment. In addition, he 
stated that more immediate gains in forecasting storm move- 
ment would come from examining the surrounding synoptic 
meelds than from detailed modeling of the inner core of the 
storm. He bases this assessment on the fact that for at 
least the foreseeable future, detailed data in the immediate 
Vieimity Of the storm will not be available. Furthermore, 
with continued advances in satellite sensors and satellite 
interpretation, data will be available to better define the 


Piraeeonment around the storm, if not the inner circulation. 
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This research is an initial study to determine the feasi- 
Belity Of applying a synoptic Map-typing scheme to the prob- 
lem of forecasting tropical storms in the Western North 
m@eliie Ocean. The basic concept is that tropical storm 
movement is directly related to the forcing of the general 
omeic height field in which it is embedded, and that this 
forcing can be represented by a set of map types. Because 
of time limitations, the techniques applied could not be 
Meeemigea. Rather, the aim is to demonstrate that the basic 
concept is sound, before perfecting the technique and adding 
the detailed refinements typical of other analog schemes. 

Analog techniques typically involve relating the current 
storm with historical storms that are temporally and spva- 
tially similar. Among the parameters used in the Atlantic 
Ocean by Hope and Neumann (1970) and in the Pacific Ocean by 
Jarrell and Somervell (1970) are distances to the subtropi- 
cal ridge and to the nearest trough. Instead of deriving 
such subjective parameters, this study directly uses the 
synoptic fields to objectively determine which relationships 
are the important ones. Tse (1966) attempted a similar 
Study, in wnich he subjectively developed five groups of map 
types. He successfully used these to stratify his sample 
cases prior to developing a pressure gradient nomogram to 


Peeadtet typhoon movement. 
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Meehis scucy, analyzed W-vaiue (deviations 
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Gard neights at specific pressure leveis) fields were ob- 
metmea from Systems and Applied Sciences Corporation, 


Ponperey, CA. A total of 504 individual cases, fro 
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aa 1 enu- 
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ary i1S67 to 31 December 1976, were intervolated to a stan- 
@eear grad (see Fig. 1). Since mav types must be constructed 
trem SyPcptic situations relative to the tropical storm, 


14 had to bewmovable in space so that the storm is 


- 


eee ys loceted im the same position. In eacne 
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megane) 5000 mb rields were utilized to determine whicn 
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Specific level) prcvided the best steering information. 
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Empirical Orthogonal Function (EOF) analysis was used 
Meemar® FSasOns. First, it can be used to separate the ran- 
Bom mOise from the true synoptic patterns. Secondly, it can 
De usec tc reduce the number of descriptors required to de- 
fees ene Synoptic field from the 120 grid points to a 
smaller number of EOF coefficients. Map types were selected 
using least squares statistics from the EOF representations 
Beene P-value tields surrounding the 504 storm cases. The 
basic procedure is analogous to that of Lund (1963). Each 
of the historical fields is then screened a second time to 
assign it to one of the selected map types. 

As a preliminary test of the usefulness of this map 
typing, a forecast position is derived using archived best- 
[meee <« slOrm poSitions. stach forecast storm position is 
derived from an unweighted average of best-track positions 
(rotated so that each case has the same movement for the 
previous 12 hours) for all the other storms within the map 
type. This procedure is applied to all three levels and 
Pee e-Orecect time intervals. The results are then used to 
Support the hypothesis that storm movement is related to 
the surrounding synoptic height field. Finally, several 
Suggestions are made for further improvements to this tech- 
nique, and for using the basic concept to develop another 


type of analog scheme. 
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iiveeeont eh ACOULSTIT ION AND EMPIRICAL 
ORTHOGONAL FUNCTION ANALYSIS 
Pee of ORM SELECTION 

Peer rom 4a ten-year period from 1 January 1967 to 31 
December 1976 form the sample. All storms between 1 October 
ico ene 31 December 1972, and 7 September 1973 to 31 Decem- 
ber 1373 had to be omitted because the D-value fields were 
Memweavaillable. Only two historical data sets, the best- 
track storm data and the upper-level D-value analysis 
fieids, were required. Both sets were available from the 
Fleet Numerical Oceanographic Center's (FNOC) archive files. 

The actual selection of the 504 case studies was made 
Smecie basis of the following criteria: 

(1) The storm had to be at least tropical storm 
Beeength, with a maximum wind speed in excess of 35 Knots; 

@oeelie storm had to have at least 12 hours of prior 
meeeciana persist for at least 30 hours; 

Gomeeonl, ©0 GMT and 12 GMT initial storm positions were 
considered since upper-level analysis fields are only avail- 
able at these times; 

Mpmeeach Case had to have ali three levels (850, 700 


SaaGmoolemp>) of data; 
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(3) To insure independent samples of gridded D-value 
maps relative to the same storm, each case had to be at 
Meect oo MOUrS apart. However, this restriction did not ex- 
eeude data trom a second storm occurring on the same date, 
Decause each case would be gridded according to its own 
location, and would therefore represent independent samples; 
and 

Sees’ storms located south of 10° N lat were elimi- 
mated since that would extend the grid into the Southern 
Hemisphere. 

A standard grid was required that could be relocated 
Deased on the position of the tropical storm. Neumann (1979) 
presents a grid that has been frequently used in Atlantic 
tropical cyclone forecast models. We selected a similar 
feels grid and positioned the storm center at the 70 posi- 
tion (see Fig. 1), since most storms have a west-northwest 
movement. The storm center and grid were located in the 
FNOC 653 x 63 hemispheric analysis fields and the D-values 
were interpolated to the grid points using the BSSLG sub- 
Beeline from the FNOC computer library. This subroutine 
interpolates a data field of dimensions (M x N) for a point 
marmoesessel's central difference formula. 

Best-track storm positions are collected every 6 hours, 
beginning 12 hours prior to the forecast time and extending 


to + 84 hours. Since upper-level data fields are only 


i 





available at 00 GMT and 12 GMT, map typing can be done only 
twice each day. Because the analysis fields are not avail- 
able until about 03 GMT and 15 GMT, a 24-hour forecast from 
map time would only be an 18-hour forecast by the time the 
Joint Typhoon Warning Center (JTWC) could use it to prepare 
a Warning. On the basis of these considerations, forecasts 
are made at 30, 36, 54, 60, 78 and 84 hours from analysis 
time, with new forecasts available for the 06 GMT and 18 GMT 
Warnings (see Fig. 2). In Table I, the 504 cases are broken 
down by the number of cases in the sample with sufficient 
best—track data to be useful in developing the analog 
forecast positions. 

An alternate approach of using the most recent upper- 
level analysis with the current storm position (+6 or +12 
hours from analysis time) was considered unsatisfactory for 
two reasons. First, there would always be a discrepancy be- 
mveoch tie present storm position and the actual position of 
the Low pressure center in the map. Elimination of this 
Giscrepancy would require smoothing in the region of the 
peVvomsstOrm DOSItiOnN and the introduction of a bogus storm 
at the present position. Secondly and most importantly, the 
Se iemeane COMdItions to the north of the storm that may affect 
storm motion are generally moving in the opposite direction 
of the tropical storm (see Fig. 3). Therefore, the synoptic 


pattern north of the storm would generally be translating 
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scheduled warnings issued by JTWC each day. 
The X marks upper level analysis time; o 
represents time that forecasts would be avail- 
able from this research, * represents the time 
warnings are issued; and $ represents warning 
forecast times. This diagram illustrates the 
reason for the 30, 36-, 54-, 60-, 78- and 8H- 
hour forecast intervals. 
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Number of cases available ver forecast time interval 


Porecast Interval Number or Cases 
30 504 
36 482 
oF 3.07 
60 360 
Jas (ome 
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normalized at each point to use the EOF method of extracting 
"signal" from "noise". This "noise" comes primarily from 
errors in the observations used to construct the D-value 
fields and from subgrid scale phenomena. The normalization, 


performed by rows, is 
= —1___J Oe 


where 1 represents the map number; j represents the grid 
Boint ; Ass represents the standardized D-value, Ks represents 
eaenmeanm Cover i) of the grid-point values j; S 5 represents 
mae wstandard deviation at grid-point j; and oe represents 
tne D-values of the ith map at grid-point }j. 

Two other normalization procedures were tried: 1) nor- 
malizing each map (i.e. column), using a and s; and 2) nor- 
malizing over the whole matrix using x and s, as constants. 
Fach method emphasized certain aspects of the data field at 
the expense of another. Fig. 4 and Fig. 5 are examples of 
the results of the different normalizing schemes after the 
Bol Comemcations and truncation. Normalization by grid point 
Seeder) insures that the variance at each grid point 
is given equal weight and emphasizes the features in the 
southern portion of the maps. Normalization by map insures 
that the variance of each map is given equal weight. Since 


most of the variance, within each map, 1s usually in the 


alge 
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northern region, this method generally emphasizes the north- 
peg@ereeions, but not always. Normalization by the entire 
geld insures that the variance of each point of each map is 
given equal weight. The effect of this method varies from 
map tO map and depends on how the variance of the map is dis- 
merpuced relative to the average variance over all cases. 
mame oP and Fig. SA;"the actual map data, rirst note 
Bie decail and then the "noise" in the grids, primarily 
represented Dy the unsmoothed isopleths. Then note how much 
smoother the truncated eigenvector fields are regardless of 
ifeemormalization scheme. In Fig. 4B and Fig. 5B, the grid- 
point normalized maps, the definition and strength of the 
lmmeougns in the westerlies in the northern portions of the 
@rid are poor. However, the tropical sections of the grid 
are represented well. Particularly note the low pressure in 
men southwest corner of Fig. 5B that was poorly analyzed in 
Pema ric. oC and Fig. SD. The maps normalizea by the entire 
map's mean and standard deviation, Fig. 4C and Fig. 5C, have 
sood representations of troughs in the westerlies. However, 
they miss the intensity of the low pressure area in the 
Beterwest Section of Fig. SC. Almost the same comments can 


~~ 


Peieeoeooer Fig. 4) and Fig. SD as for Fig. 4C and Fig. $C, 


On 


The method of normalization by grid points used here (eqn. 
(2.1)) was selected to keep the large variance in the north- 


Penmperr of the grid from overpowering the variance in the 


oe 








eeeemern part. Atso by following this method, follow-on 
studies will be able to use the Preisendorfer and Barnett 
mee 7) method for truncating the EOF series. In Fig. 6, 
Bie maps of the 850, 700 and 500 mb gridc-point means anc 
standardc deviations, that are used to normalize the grid 
fea (2,2), are presented. Particularly note the north- 
Beuel variations of the means and the large standard deviae- 
@eons in the northern portions of the maps at all levels. 
Mrere 1s also a relative maximum in the standard deviations 
mee me Location of “the tropical cyclone, which indicates a 


eaowe Variation in storm intensity from map +o map. 
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See er LRICA:i, ORTROGONAL FUNCTION (CEOF) ANALYSIS 
The concept of using EOF analysis in meteorology was 

mest introduced by Lorenz (1956). His deseription formed 
[ee Desis For a wide range of uses of the technique in the 
MeecoOrological fieid. Kutzbach (1967) provides a concise 
Meaenematical and descriptive outline or the technique upon 
Sey nis research is based. The object of using the EOF 
mo oaneaguie 25 twO-rold. First, 2 is used in an attempt to 
u ‘noise" in the D- 


nal" from random 


co 
— 


f4- 


s 
ields. Secondly, it provides a method to recuce the 
number of variables from a large number of highly correlated 
variables to a smaller number of independent variables. 


Only a general description of the concept and procedure is 
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The D-value means Cleft side) and standard 
deviations (7teh: sade es ror 50 ma (hottcm).. 
700 mb (middle) and 500 mb (top). The 
isopleths are D-values (m). 
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included here. A more detailed mathematical description of 
the process is presented in Appendix A. A listing of the 
computer subroutines used in the analysis is given in 
Appendix B. 

A 120 x 504 matrix X was formed, such that each grid 
point in each map is normalized in relation to that grid 
point over the 504 maps. Each map may be thought of as a 
vector in 120-dimensional space. The 504 vectors represent- 


ing each map will tend to be highly correlated. Therefore, 


}-- 


1t 1s expected that the points representing the position of 
the head of the 504 vectors will tend to cluster in space. 
Miemebjeetive of the EOF analysis is to use the 504 points 
EO determine aeset “Oc l720serenogonal vectors to represent 
the basis of the 120-dimensional orthogonal space. To de- 
termine this 120-dimensional space, first form the symme- 
tric normalized covariance matrix. The columns are vectors 
representing the variance in the data matrix. The first 
aXis in this new space, hereafter referred to as the first 
elgenvector, will be determined such that it will be a li- 
near best fit to the variance vectors. This procedure is 
equivalent to maximizing the explained variance (eigenvalue) 
and minimizing the unexplained variance (residual) of the 
data by the first eigenvector. The second eigenvector will 


aeeeaeecad Linear best fit to the variance vectors, with the 


additional requirement that it be orthogonal to the first. 


as) 








mms, chis eigenvector will explain as much of the remaining 
Meeeiamee as possible. The third eigenvector will again be 
meee linear best fit to the variance vector, but it must now 
be orthogonal to the first two eigenvectors. The procedure 
Will be continued until all the variance is explained, which 
will be equivalent to having an eigenvector for each of the 
m0 dimensions. 

One may think of EOF analysis as being analogous to a 
Fourier series. Instead of representing the data as a linear 
combination of space-orthogonal sines and cosines, the space- 
orthogonal functions in the EOF analysis are derived from the 
actual variance in the data. As with Fourier analysis, each 
eigenvector has an associated (time) coefficient. This co- 
efficient, when combined with its eigenfunction and summed 


over the entire eigenspace, provides an exact reproduction 


of the original field. That is, 
M 
DG a 3p) ee » } EOF (x) Cnt t?) CZs 
m=1 


where X represents the value at grid-point x of map t in the 
data matr1x; EOF. represents the empirical orthogonal func- 
tions which are derived by relating the unit eigenvectors to 


mn 


the original 120-dimensional space; and ©. are the time 
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coefficients (in this case time is determined by the differ- 
ent map times) related to the m th eigenvalue, and M = 120. 
Miese cOectficients are time orthogonal (Kutzbach, 1967) and 
are obtained from the eigenvectors and the original data 
field. The specific mathematical relationships are described 
in Appendix A. When a linear combination of elgenvectors 
Bim@e their respective time coefficients are translated into 
the original space, the eigenvectors represent a series of 
functions in the same manner as sines and cosines represent 
Mmer Loins in Fourier analysis. A truncated linear combina- 
meiotic. (2.2) with M less than 120) ean be shown to 
approximate the original maps with the "noise" removed 

Bee vsendowrer and Barnett, 1977). 

Upon completion of the EOF analysis, the eigenvectors 
and eigenvalues represent solutions to the maximization of 
Mmecermlained Variance" problem. Thus, depending on the 
geophysical parameter being analyzed and the quality of the 
data available, the first few eigenvectors often represent a 
large fraction of the variance in the data set. Preisen- 
dorfer and Barnett (1977) have developed an empirical 
Hommila for separation of true "signal" from "noise" in EOF 
analysis. Since each eigenvector represents a weighted por- 
tion of the variance, this separation can be accomplished 
by truncating the number of eigenvectors included in the EOF. 
Because of time constraints for this study, an arbitrary 


moii@amlon was Made at eigenvector 10. Without formal 
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application of the Preisendorfer-Barnett techniques, we can 
not be assured of the statistical significance of the higher 
Order eigenvalues. This is one aspect that needs to be 
optimized in any subsequent studies. 

Because most of the variance is contained in the first 
rew eigenvectors, excellent "noise-free" representations can 
be made from just a few of them. The cumulative percent of 
the variance expressed by the linear combination of eigen- 
vectors derived from the 504 case studies is shown in Table 
Ti, For example, 83-85 percent of the variability is ex- 
pressed in the first 10 eigenvectors at all three levels. 
Maeeact, beyond a certain point, the eigenvectors represent 
such a small portion of the variance that they just repre- 
Sent noise in the map fields. In this manner we can reduce 
tne number of variables required to represent the maps from 
meeeerid points to 10 coefficients of the orthogonal eigen- 
Meetors. In the next section, we use these 10 coefficients 
eect) ) to derive a series of characteristic maps which we 
will call map types. 

The maps in Fig. 7 represent case number 1 and are an ex- 
ample of how the linear combinations of the EOF build the 
map. The last map in Fig. 7 is the same as shown in Fig. 4A. 
Withmout describing in detail each map in Fig. 7, it can be 
seen that the more drastic changes in the pattern occur in 


—MmOeMenocti of westerlies after the anclusion of more than 10 
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dives) oes JEL 


Cumulative percent of variance and the actual variance 
oe.) OP elgenvalue as expressed by the linear combina- 
tion of eigenvectors for each level. An arbitrary 
cutoff at 10 eigenvectors was used here. 


Elgen 850 mb 700 mb SO mers 
vector Percent Variance Percent Variance Percent Variance 
1 34 40.15 tis PaaS ao) 35 39,36 
2 o1 Zl S 46 Dono ak Haale sal 
3 50 A es Oo 6 lee iey 59 Slee) 
+ Sis 8.42 eZ peo 65 bats 
S a2 6.09 68 Te colli 70 ba Us 
6 aS 4.14 72 5.14 75 Sao 
tl 78 S20 42 76 See 78 ey hs 
8 81 een Sle. 719 Sie) S, 8 1 Coal 
g 8 3 22756 ol 2 ow 83 2.43 
Ke 85 papas 83 Pes) 85 Peealealt 
oi). 86 sire el 84 1.98 86 Ie se 
lee 87 reo 86 iesc a 88 Le 
3S 8 8 eats 8 8 dei te A, oo ieee 
14 che) 20S 89 sie e 20 hl 
ilo 90 Sas go ecole on ie OS 
ike ik 3S Sh IB Bis SAL SiS 
17 a 84 Sul 97 a2 78 
18 oe ee a2 74 os a 
ine Shs Su a2 62 SIs 69 
20 93 55 2 61 Q4 58 
m7) 0.0 Eo 0 co 
tc 15 terms. With just one term, the basic flow is estab- 
lished - moderate westerlies, weak gradient in the tropics 


and a trough where the storm is located. After including about 
15 terms only slight refinements of the pattern take place. 
The "noise" in the map starts to become more obvious with the 
inclusion of 40 terms. The isopleths begin to show irregu- 
larities and are not as smooth as in the earlier maps. Be- 


yond 80 terms it is almost impossible to see any differences. 
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Figure 7. An example of the effect of truncation of the 
POrerowwoul sou , 18!)Man 1956/7. - Tae number above 
the map is the number of terms included (the 
truncation point). The isopleths are D-values (m). 
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IIL. MAP TYPING 


Bemalris point, three choices are available. The EOF 
peatysis could be treated just as an approximate method to 
Seemenate the noise from the data. Then "noise-free" maps 
Berea be produced and map typed exactly as in Lund (1963). 
Mae Sseconc choice, and the one followed here, was to use the 
truncated EOF analysis to compare the maps in eigenvector 
eae. ihe actual map typing procedure is similar to the 
Me DmOcCeduUre, except that instead of having 120 variabies, 
enere are only 10 in the truncated eigenvector space. 

Peemira possibility is to consider the eigenvectors as 
the map types. The grid-point normalized fields could be 


Begeared in either least squares approach or a direct cor- 


f)) 


Me 8er2On approach, to the actual eigenvectors. This method 


Fh 


Meme Gerenosen in this study Since the statistical signifi- 


cance of each eigenvector was not known (Preisendorfer and 
Barnett, 1977). By using the coefficients, and therefore 


essentiaily the approximated data sets, we hoped to minimize 
the problem of not knowing the statistically derived 
matmecdr: 107 DOint. 

Pemeecseme pea orevicously, the eigenvectors represent tne 
Beeeminecpace. tt is the coefficients that actually locete 


Mee DOsdulon Of each individual Map in that space. Although 
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each map is an independent sample, the maps are highly cor- 
related and will when plotted in space tend to form clusters. 
The idea is to find the map that is most representative of 
poe cluster and call that the map type. To do this, each 
map 1s compared in the least squares sense to all the other 
Maps. This is done by using the coefficients of the trun- 


Seod Gigenvector space in the following manner 


JEG 


: > ie 


1=1 


where 22 ay, is the least squares value comparing map j with 
map k; ie is the coefficient of the ith eigenvector of map 
ae Cok is the coefficient of the ith eigenvector of map k: 
and pe is the variance accounted for by the ith eigenvector 
[Mice truncated eigenvector set. The ie value for each 
map with all cther maps is computed and stored. 

Now, search for the map that 1s most representative of 
mes kientest” cluster and designate this as map type 1. To 
do this, the map that had the largest number of other maps 
Falling within the least squares radius of five units was 
designated map type 1. The least squares radius of five 
units was selected arbitrarily on the basis of a trial and 
error method. Since map type 1 was the most common map, it 


should have the largest number of maps associated with it, 


but no so many that the sample is so depleted that the 
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other map types can not be determined. The least squares 
radius of five units appeared to give the most representative 
sets. It 1s possible that a more rigorous approach for find- 
ing a statistically valid value of the cluster radius could 
be determined. 

Type 1 and all maps selected within the least squares 
radius of five were removed from the sample prior to the 
oeerton of the next map type. Since it is anticipated that 
mieemext Cluster would not be as tight as the first, the 
least squares radius was increased by two units. Map tyve 2 
is then selected as the map with the largest number of remain- 
ing maps within a least squares radius of seven units. Then 
map type 2 and all maps within seven units of it were removed 
from the sample. The process, with an increasingly larger 
radius, is continued until at least 470 maps were selected. 

Finally, all the remaining maps are lumped together in 
the last type. The actual map type is selected in the same 
manner as all the others, except that now all the remaining 
maps are placed with this type regardless of the value of the 
least squares radius. 

When this process was completed, there were 11 or 12 map 
types for each level, each with a number of other maps simi- 
mieten it. AS the typing process evolved, there was no 
attempt to keep the least square radius circles from over- 
Zepiceweenus, che Map typing is not complete, because a 


Map in type 1 may actually be closer to map type 3 than it 
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is to map type 1. Therefore, each map must again be compared, 
in the least squares sense, to each of the selected, most 
representative map types. Then each map is grouped with the 
mye it 15 nearest to in space. Now the map typing is com- 
piete. Only the results of the 500 mb map typing are dis- 
GCussed here. Similar comments could be made for the 850 and 
POG mb map types, but in the interest of brevity just the map 
types and forecast results are presented in Appendix C. 

Maem Or the eleven 500 mb map types in Fig. 8 is clearly 
Meee pcnt. The following is a discussion of the distinguish-— 
ing synoptic features for each map type: 

Map Type 1: The storm is situated moderately deep in the 
meopitcs. fo the west of the Serum! them pressure sarad lent is 
very weak. The broad ridge to the east of the storm is also 
weak. Well to the north of the storm the pressure gradient 
of the westerlies is weak PRES erouenent indicated to the 
north-northwest of the storm at a distance of approximately 
foo n.mil. 

Meapeelyoe 2: Here the storm is deep in the tropical belt 
with a strong east-west subtropical ridge centered aporoxi- 
a weeoU0nm.mi.L. almost due north of the storm. A weak ex- 
tension of this ridge separates the storm from the broad low 
molm@nem. to the west of the storm. To the east of the 
storm the influence of the strong subtropical ridge extends 


tO the southern boundary of the map. 
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melee 6, 





The 500 mb map types. The number above 
each map refers to the map type number, 


1.e. in the order selected. Isopleths 
are D-values (m). 











Map Type 3: In this case the storm is near the north- 
ern boundary of the tropical region, with weak pressure 
gradients to the west and west-southwest. The center of a 
strong cell at the eastern extremity of the subtropical 
ridge is evident approximately 1200 n.mi. to the east-north- 
east or the storm. The strong westerlies that begin only 
600-700 n.mi. to the north of the storm further intensity 
to the north of this ridge cell. There is a weak trough to 
the west-northwest. 

Map Type 4: Here the storm is only 300-400 n.mi. from 
mmemcotrong westerlies. It is located between two cells of 
P@emouberopical ridge with the stronger cell due east at the 
edge of the grid. There is a trough in the westerlies almost 
due north with the storm apparently at the southern extent 
Seeyche trough. 

Map Type 5: Again the storm appears to be located deep 
in the tropics because of the surrounding weak pressure 
gradient and the absence of the strong westerly flow in the 
northern portions of the grid. There is a broad weak ridge 
memrme east of the storm. An extension of this ridge can be 
seen in the westerlies far to the north of the storm. 

Map Type 6: In this case the presence of the storm is 
barely discernable. The implication is that the storm is 
elther very weak or that its influence is not seen at 500 mb. 
Another point that needs to be made here is that the data 


on the standard grids is interpolated directly from the 
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analysis fields as they appear on the FNOC data tape. Some 
fields may have the storm location and intensity inserted 
manually while others may not. In the scope of this re- 
search, there is no way of knowing whether or not the bogus 
Meee were applied to any particular map used in our data 
set. Ihe center of the subtropical ridge is positioned only 
about 500 n.mi. to the east of the storm with its influence 
extending well to the west of the storm. The westerlies to 
the north are of only moderate intensity with a trough-ridge 
Pageerm in the northwest quarter of the grid. 

Map Type 7: Again the storm circulation is only weakly 
apparent at 500 mb by the slight deformation in the south- 
west section cf a strong subtropical ridge to the northeast 
of the storm. In this case however, there is a strong mid- 
Mretrude trough system extending to the latitude of the 
storm to the storm's west. The westerlies through the 
trough are moderate, increasing in speed nortn of the sub- 
Mmeopical ridge position. 

iememope SG: In this case the storm's influence is felt 
strongly at 500 mb. The D-value at the storm's location is 
less than 240 m, which is more intense than all other map 
types except map type 1l. The storm is located between two 
Besememeell= in the subtropical ridge. The pressure gradient 
to the east of the storm is much steeper than to the west 
within 600 n.mi. of the storm. The westerlies appear to he 
of moderate intensity with a well defined trough to the 


meoren=—northwest. 
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Map Type 9: All indications here are that the storm, as 
represented by the trough in the westerlies,has almost en- 
tirely departed the tropics. In the western portion of the 
grid the subtropical ridge has either completely broken down 
Seeene Storm is so far north that the western part of the 
mage 15 south of the grid. The ridge center to the east of 
the storm is strong, and there is a strong westerly current 
Bemeetes Orch of it. There is a trough in the westerlies 
eemecst Cue north of the storm. 

Map Type 10: As with map types 6 and 7, the storm here 
1s barely discernable at 500 mb. It appears here as a break 
in the ridge. The ridge axis is oriented from east-northeast 
©f the storm to almost due west of the storm and is much 
stronger in the northeastern extension. The westerlies are 
strong with a trough located north-northwest of the storm. 

Map Type ll: In this case the storm appears deeper than 
in all the others. The storm appears to be 700-800 n.mi. 
pemieneor ene subtropical ridge axis. A strong cell in the 
ridge is evident to the northwest of the storm. A mid-lati- 
tude trough is evident in the northeast corner of the grid, 
and a trough extends through a weakness in the ridge to the 
storm. The westerly flow 1s well removed from the storm 
Ole o-=1700 n.mi.). 

Pitheonis briet Gdtseussion Of Gach Map, If is clear that 
each represents a different synoptic situation. Fig. 9 and 


bua. 10 show examples of map types 1 and 2 storm cases from 
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Race 9. 





Examples of map type 2 (500 mb). The letters 
above the maps are for reference for the text. 
The isopleths are D-values (m). 


ul 








Figure 10. Examples of map type 1 (500 mb). The letters 
above the maps are for reference for the text. 
The 1sopleths are D-values (m). 
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the sample (more examples of the other map types are shown 
eepoendix D). Sinee each map had to be classified with 
ene of the characteristic maps, the least squared distances 
between the individual maps and its characteristic map may 
Vary considerably even within one map type. Thus, it is to 
be expected that some maps show greater similarities to 
Beeiny Characteristic map than others (see Fig. 8). 

In Fig. 9, the map type 2 maps have a strong ridge to the 
Meren of the storm. The actual position of the primary cell 
varies from north-northwest of the storm to northeast of the 
eeoom. All the maps show the storm deep in the tropics. 

The intensity of the storm varies considerably. In Fig. 39, 

A and E, the storm is very weak at 500 mb, while it is very 
Seems in Fig. 9, B and D. The westerlies, well to the north, 
vary in intensity but are generally fairly weak and the 
pattern 1s similar to the map type. All the cases indicate 
a broad area of low pressure in the tropics due west of the 
storm, except in Fig. 9C when the low pressure is to the 
west-southwest. Therefore the basic patterns are very 
semilar. | 

In Fig. 10, the map type 1 examples do not show as much 
Similarity to their map type case as do the map type 2 ex- 
amples, but the basic patterns are still evident. All maps 
show strong evidence of the trough in the westerlies to the 
north-northwest of the storm that is found in the character- 


istic map. All cases show westerlies that are weak to 
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moderate in intensity. The storm is always moderately deep 
iiente Eropics with the exception of Fig. 10E. In oe peely Ol. 
B, C, D, and E, the storm is weakly represented at 500 mb as 
meme Map type storm. All show the ridge basically to the 
PeernOr eRe Storm. Again, it is not too difficult to see 
the similarities. 

Arter examining Fig. 7 and noting that map type number 1 
may have needed more terms to adequately define the details 
Or the synoptic features, one actual size (500 mb) with the 
EOF series truncated at 20 terms (approximately 93 percent 
of the variance) was tried. The map types selected using the 
20 EOF terms are shown in Fig. 11. The selection procedure 
was Similar to that described above except that a larger 
(15 versus 5) initial least acjueiaes distance was required 
initially in determining the characteristic map types. The 
inclusion of more terms appears to better separate the maps 
in elgenvector space. The maps have considerably more detail 
and the D-values in the troughs andridges much more closely 
resemble the actual maps. It is clear that the map typing 
is very sensitive to the number of EOF modes that are re- 
Tained in the analysis. It 1s possible that 10 EOF terms 
may not be sufficient to properly represent the significant 
features for several possible reasons. First the data should 
have already been somewhat filtered, since the FNOC fields 


have been objectively analyzed. Because the grid size 
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Figure ll. The 500 mb map types using 20 terms of the EOF. 
The number above the maps refer to the map type 
number. The isopleths are D-values (m). 
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used is not that much different from the tropical grid 
size of the analysis fields, there would appear to be iittle 
subgrid scale features to be removed from the grid-point 


values. 
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IV. ANALOG TRACK FORECASTS BY MAP TYPES 


With each map now typed, a forecast position is derived. 
The following procedure is simple and later can be greatly 
refined, as indicated in the next chapter. As indicated in 
the introduction, the object here was not to optimize the 
Pecimteue, but rather to see if a representative forecast 
could be made based on the map types. 

Betore any Comparison of forecast positions could be 
made, all the positions had to be standardized. The O-hour 
@estelon, Or map time position, of the storm is converted to 
Mmieeiand I40° FE. This position 1s arbitrary and does not 
atfect the results. It only serves as a reference position. 
All other best-track positions are then converted to a lati- 
tude/longitude position relative to this reference position. 
Then, each case is selected in succession. All other cases 
in its same map type serve as forecast analogs. Before the 
forecast positions are actually derived, all the best-track 
positions are rotated about the O-hour storm position so 
that each analog has the same direction of movement for the 
12 hours prior to map time. There are no restrictions placed 
@Methis rotation. The equations that accomplish this 
motamileon are , 


CALA SLALT + (CRA —-PSmatj es FCSsT/12 
CAcON = SLON + CPA EON-PslON) <« Fest/i2 
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where CALAT and CALON represent the rotated analog latitude 
and longitude; SLAT and SLON represent the unrotated analog 
latitude and longitude; PALAT and PALON represent the prior 
Be-nour analog latitude and longitude; PSLAT and PSLON 
represent the prior 12-hour current case latitude and longi- 
tude; and FCST represents the forecast interval, i.e. 30, 
36, 34, 60, 78 and 84 hours. It also assumes that changes 
in latitude and longitude by an equal increment represent a 
meue distance. 

All rotated analog storm positions within the same map 
type (excluding the map for the case being forecast) are 
then averaged, and this average storm position becomes the 
forecast Dosition for that case. These forecast positions 
are compared against the actual standardized best-track 
position. A mean error and standard deviation are computed 
for each map type and for the full 504 case sample. A sum- 
mary of the results for the full sample is provided in 
Table III. See Appendix C for the map type summaries for 
moO and 850 mb levels. 

It was hoped that this study would show that the easy to 
Mmereecast Situations Could be separated from the difficult 
Ones, and a quantitative estimate of forecast accuracy could 
be made. In a limited sense this can be shown in Table IV 
Pyecombparing Map type 2 and 3. The mean error and standard 


deviation for map type 2 are less than those for map type 3 
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summary of mean forecast errors (n.mi.) and 
standard deviations (n.mi.) for the analog- 
type forecast. 


Fest 850 mb 700 mb 200 mb 

bg ws Mean Std Dev Mean Std Dev Mean Std Dev 
30 S67 aes gas 137.9 Sse 2 eo OZ Shome0. 

36 Mo. 11g .7 A Se. 2. 6 Oe be ORO ae eS 

54 Boe Ss OLB. 1 Deere poco: Gece! 26474 182.5 

60 Sore. 20742 OFS Suet BZ eS ak Gaal ove tay. wae sed re: 

78 woo.4 6287.2 Oe eee 2 8 7a SZ Ue 220859 

84 fee -317..:5 OSes 321. see Oa 

im every forecast time interval. This would indicate that 


storms in map type 2 are easier to forecast (by this analog 
procedure) than those of map type 3. Based on the maps in 
Fig. 8, most forecasters would probably agree with this 
assessment. However, a test with independent data after 
optimizing the technique will give a much better indication 
Seepchie Significance of this concept. It does show that the 
Meerential exists. 

It is obvious from the average errors, as presented in 
Beeeeseiil and IV, that the results do Not represent any 
major breakthrough. However, it must be remembered that the 
techniques have not been optimized and many refinements (see 
Chapter V) can be applied that should give better forecasts. 
In any case, the errors are not that bad when compared to 
the official forecast errors and represent a solid base from 
which to build. The official forecast errors for the speci- 


Eeewesses in this study are not known. A rough estimate 
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of the official JTWC warning error was computed by taking 
the mean error for each of the years included in this study, 
summing them and dividing by the number of years (10). The 
results indicated an average error of 112 n.mi. at 24 hours, 
aooeae mi. at HOES Inet enema tv nouns. These 
ewemage errors include all warnings issued for storms of 
Pmaompreal Storm intensity for the time frame involved. 

Another question regards the level that gives the best 
aeeeming motion. It is obvious from Table TIL that the 
data from all three levels produced remarkably similar 
errors. Why this occurred is uncertain, but it is con- 
jectured that the adjustment made for the past 12-hour 
motion strongly influenced these average error statistics. 
To further examine how the schemes at each level performed, 
a block of 20 storms was selected at random from the data 
sets and presented in Table V. The data at each level for 
these storms is remarkably consistent, which may again be 
the result of the corrections made for the past 12-hour 


movement. 


nN) 





eae l/r mt 
eer) OL 
Smegedh, oD 
Om GO 
Be Leek7 0 ae. 
Gage de S “eal 
Cee es cl 
hee Sy iam oa | 
Oily res) 
CO br 
Soles 1 eh all 
foe 7 O aL 
ee 7 baer 
ee nino", SL 
Stir Saab 
e iil, 0 el 
oe ot ob 

Cee le 8 ab 
Gece /-S. Unt 
Gailey | lige 


TVALOV 


SU a leOcd iii y, 


LO LO LON LON LON ON ON ON OLN 
AON NN MEN ak AA TN 
Nee re Nee Nee ee NY Ne Ne OS 


“—_- 
Pcp | | 
Oo So 
Ww aw 


Ct) 
Gi) 
(2) 
Ce) 
Ce) 


Seoc ly « irl 
Celle yee 
OC a si 
Se lhy os 0h 
Ge Sel. t 
Oe CL Ge 
Creo Gilead. OE 
SPO cally eee sib 
Cy ott 
0-6 e/a | 
e Ge / cel 
Cee eal 
lye tole sé asaya 
tie Siew o. aa L 
OF Cm) 3 hit 
eee mm/h ll 
£ Ginu/<c 6 

Gasiity de cr 
Oe ec ly oll 
I> Cin, iele 


{WN OOS 


"qu 00S PUP OOL 


Fa a a a 
WHwmwmwonm 
Nee Nee ee NL Oe” 


oe 
Gist) 
— CT) 
C1: 
ly) 
Ct) 
Geo 
Ct) 
(ti) 
(2) 
(1k) 
Ce) 
(2) 
(ti) 


eel) Co eree IL 
SOB as IIE 
OF Ghee 
BCG cally neal 
Coca Oe c LE 
Sek ae el 
C79 el 0 
cr qe ey eal 
FeO Aly Sle 5 IE 
Ob og 01 
Gb ely tc 
OP ES i/o, aed 
STi os Sully ee pea SyaE 
See, Giirc il 
Ce Clue Cork 
G00 ly ec. 
Gila 6 

Seas Es ac ae a ke 
‘ieee onl Bey 4m al 
Oo -CWn/ ie cl 


dW O04 


~o- 

CN ON OE OD OR ON OE LR LI OI OE OTS 
HMO DA TA TA MTN NN ON N 
Nee Ne Nee ee Nee Nee “Ne NL NL ON Se” 


o™ 
ae 
ast - 
wt 


Co 
Ce) 
Gre) 
eo 
Coy) 


Cab: Ly earl 
Ga tly aD 
oc ei 
see e Ieee Dull 
S) te llyers. 764 
US scl Ome 
eee Cae al 
Gi Ie leah a 
(om OMe elk 
Gee E700 
Weigel si GIL 
8) eedby (3 tei 
Sts ee pal gS) 
67 St 1/8 al 
eu ile 
Bred by) (9 el 
Sea Ur a6 

Tees alee EL 
C2 Ee 0a 
ee Ly 3h ane 


dW O0S8 


hoE 
COE 
C6E 
Jee 
O6E 
Bee 
88E 
Eee 
Oe 
Sits, te 
hge 
etere 
CBE 
Pee 
ORES 
BEL 
Cie 
Le 
Jie 
neo 


YdeWAN ASVO 


*(dequnu edAQZ dew) SuoT/jAeT se usats 
pue (7 oOtTt ‘N oO0T) UoTZTSod Unoy-Q peztTpuepuej,s sy} OF BATRETSed sue 
‘oss 37e apew szsSeosdojy unoy-g¢e jo satTdwexy 


"A ATV 


24 








V. SUGGESTED IMPROVEMENTS 


Pec Gnis was am initial study to verify that the tech- 
niques would work, several procedures were not optimized. 
The first procedure that needs this optimization is the EOF 
truncation procedure, using the Preisendorfer and Barnett 
me) method. In Chapter III we noted that the 20 term 
EOF map types show considerably more detail than the ones 
attempted with 10 terms. Whether this additional detail is 
Sociale Or mot (i.e., whether it is "signal" or "noise") 
remains to be proved. We suspect that the motion of the 
tropical storm is primarily dependent on the large-scale 
synoptic situation, and the 10 term EOF map types may have 
represented these large-scale features well enough to dis- 
tinguish the map types. Unfortunately, it is also possible 
mae some Of the first 10 terms of the EOF may actually be 
"noise". Thus, there may be too many terms in the EOF. 
Table VI contains the 30-hour forecast results using the 
20 terms. The average error for the entire sample does not 
show any improvement, but it does show a wider range of 
Beverage errors for the specific map types. Therefore, by 
using more terms it may be possible to better separate the 
easy forecasts from the difficult ones. Remember that this 


suggestion was made on the basis of the forecasts for one 


Se 





level and one forecast interval. More data must be gathered 


before this idea can be really supported or refuted. 


GABLE Vi 


30-hour forecast mean errors (n.mi.) 
and standard deviations (n.mi.) at 
S00 mb aiising 20 termssof the EOF. 


Map standard Number 
Type Mean Deviation of Cases 

al 1 2 PO Crr 1s 

eZ Niecy Looaet Ug 

c Jee) sal IEC) cea eS 

m ILS AS 68.9 30 

5 ings 0 65.9 37 

6 lab (0 OS 24 

| ers: 4 TSIEN S: IES 

8 Le Qa! 83.4 a? 

9 Us ree) MO SiS 17 

10 184.3 dOrer 4 15 

Jia Lia CNG 8 Li 

ae 236 0 eos uy 

JES (Esmee Doe ie 

14 a ec yeay Escenas 6 

1S) eS feel Tad 7 

16 ee oo 48.4 6 

any D2 Dee Ty aS 6 

18 eee eee 2 

Total le 2 ale OF «1: 504 


Also, as indicated at the beginning of Chapter III, 
Pmece ene proper truncation point has been determined an 
attempt should be made to use the eigenvectors as the map 
mpes. “(his use of the elgenvectors should work as well or 
better than the characteristic map type approach and has 
the advantage of much fewer computations. However, this 
procedure would work only for the map-typing scheme and 
would not work for the strictly analog scheme suggested 


fiers in this chapter. 
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Next, if the map-typing concept is to be Followed, a 
weighted average should be used to compute the forecast 
positions. The weighting function should be inversely pro- 
portional to the least squares distance from each map in 
the map type to the case being forecast. This procedure 
pomld=241ve more weight to the positions from the maps that 
most closely resemble the present case. In other words, 
the more similar the maps are, the more likely the synoptic 
forcing is similar. Therefore there is a higher probability 
miec tie future movements will be similar. Thus the fore- 
cast should be improved. 

The past 12-hour movements should also be compared be- 
fore using all the cases in the map type to compute a 
forecast movement. Any case within the map type that does 
not have a similar past movement, within say plus or minus 
30-40 degrees, should not be included when computing the 
Memeceaso position. The justification for this procedure is 
that 1f the past 12-hour movement is so much different than 
the current case, the forcing mechanisms for the two cases 
are most likely different. This could be due to rapidly 
changing patterns or to small-scale features in the synoptic 
Situation that the map-typing scheme did not resolve. Later 
Versions of the analog forecast scheme of Jarrell and 
Somervell (1970) use a similar restriction on permissible 


past-motion vectors in the analog scheme. 
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Another improvement to the methodology that should im- 
peew>] the selection of the map types is to devise a more 
Secaeastically significant means of selecting the least 
Squares distance than the trial and error method used here. 

Another more radical change to the procedure would offer 
fae epOcential for still greater improvements. This would 
memeomvary the size of the domain based on the forecast 
Meer val. fhus the domain would be smaller for the 30-hour 
forecast than the 84-hour forecast. In this study only one 
grid domain was used for all forecast intervals. By vary- 
ing the domain size we can select the features that affect 
the movement in the interval. For example, with the domain 
used in this study, the variance expressed by features at 
the extreme edge of the domain (2400 n.mi. away) may have 
no influence on the 30-hour movement, but may be very signi- 
Ficant in the longer time frames. The proposal here would 
be to shrink the domain for the 30-hour and 36-hour fore- 
casts, in order for the EOF procedure to maximize the 
variance in the area that 1s most likely to affect the 
storm's movement in that time interval. 

Another possibility is to revert to a strictly analog 
technique, in which each of the 504 cases is treated as a 
potential analog. Instead of having predetermined map 
mypes, cimply select the 10 best fits, in the least squares 


sense, to the present map. Then use these 10 best fits to 
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compute the forecast position. This would perhaps give a 
Detter forecast position, but it would not have the advan- 
tage of having a predetermined estimate of the accuracy of 


the forecast based on the map types. 
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Vi SUMMARY 


Pi eeyective or this imitial study was te usesthe hy— 
pothesis that tropical storm motion is directly related to 
Pew Orcing Of the Synoptic height fields, and to construct 
foemeypes 2anpound the Eropical storm that represent this 
forcing. After this forcing, represented by the map types, 
is determined, we then use it to derive forecast storm 
Pesitions. 

First, the FNOC analysis fields are interpolated to a 
standard grid based on the tropical storm position. Then 
che Maps are organized into a matrix and an EOF analysis is 
Gompleted. The EOF are then truncated in an effort to 
eliminate noise and the subgrid scale phenomena. The time 
coefficients of the truncated EOF are used to compare the 
maps in eigenvector space. A least squares procedure is 
used for these comparisons. The maps are highly correlated 
and will tend to form in clusters in the eigenvector space. 
By identifying these clusters and picking the maps that are 
Mesgerepresentative of each of these clusters, the map 
moines accomplished. The map types are then used to 
divide the 504 individual cases into groups. The maps in 
each group represent similar synoptic situations, which can 


be used to derive forecast storm positions. The forecasts 
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are developed by using the known best-track positions of 
fe seOrms in each map type and rotating them so that they 
have the same past 12-hour movement as the case currently 
being forecast for. An average position from the map type 
1s then used as the forecast. 

This research shows that the concept is valid. The 
procedures involved need to be optimized and refined. Map 
types can indeed be formed using the EOF analysis and least 
squares statistics. The forecast movement errors arrived 
at in this study, do not appear to be as good as the offi- 
cial forecast, but since the technique was not optimized, 
they were not expected to be as good. The fact that the 
errors in this research were as close to the official fore- 
cast errors as they were, indicates that further refinement 
of the technique is warranted. By optimizing the procedures 
and employing some of the refinements from Chapter V, this 
technique should be able to equal or better the official 


mOreCast. 
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RE PENDS 


Meinods for representing arrays of environmental fields 
by empirical orthogonal functions have been formulated in a 
number of different ways (for example Lorenz (1956), Frei- 
Bemeer ana Grenander (1965), Kutzbach (1967), Essenwanger 
meee), and Little (1980)). Most of the information pre- 
sented in this appendix is based on the work of Kutzbach 
(1967) as presented by Little (1980). 

Let A be an M component vector that represents the nth 
(n=1,...,N) observation. Let A be an M by N matrix, where 
the nth column is the observation vector AL: The goal is 
to determine a vector e which has the highest degree of simi- 
Mairicy to all the AL VeECrore.  MUEZ bach (Ide j explains tna 
the similarity is measured with the squared normalized inner 
product of the observation vectors AL and the vector e. 


After averaging over ali of A, the problem is to maximize 


Dies: Ce'A)Ce'A)' e'AA'te 
! ! = = Le 
Ce'A™)N /e'e =—heaus “Nele ( i) 
where the prime denotes the transpose. This is equivalent 
mommMccamiZing, 
e'Be 





where B represents the normalized symmetric (M x M) covari- 


ance matrix, and is subject to the conditions 


e'e =] Clea 


and 


Bese N7+CAtaA] ies.) 


Piarniermnore, iiezbach (1967) points out that the maxi- 


mization of (1-1) subject to (1-2) leads to the equation 
Be = eA (1-4) 


The vector e and the parameter A are recognized as an eigen- 
vector and associated eigenvalue of B. It can be shown 
moceeeseme € are orthogonal and that the hs are real and posi- 


tive. Thus for the whole matrix, (1-4) can be written as 
BE = 8 Ci 5.) 


merce 5 is the Mx M eigenvector matrix and Lis the M x M 


diagonal eigenvalue matrix. Also note that 
Eee =) al (1-6) 


It is now assumed that the elements of L and the respective 


columns of E have been arranged such that 4, is the largest 


aE 
A and }, is the next largest XA, ete. From (1-6) the trans- 
pose of an orthogonal matrix is equal to its inverse. Thus 


combining (1-3) and (1-5), the result is 


E'AA' = LN 
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Now define C = E'A, where C is an M x N matrix. It follows 


that 


Ae > Cen ee ee Clad) 
1=1 
Thus from (1-7), the observation vector a can be expressed 
as a linear combination of the m eigenvectors and n C vectors 
SemecOoer 7 1cient vectors. Ef the entire E and C matrices 
emeenused the exact data matrix A can be reproduced. 

In most applications to meteorological fields, where the 
Mate vectors A are highly correlated, it is found that a 
large portion of the variance (i) can be accounted for by 
Betaining only the first few of the largest elgenvalues (i) 
and eigenvectors (e). Thus by retaining only the first few 
terms of (1-7), an excellent representation of the observa- 
tion can be obtained, with the rest of the terms being con- 
sidered noise. Thus if an a” is defined as the approxima- 


jiiem of the nth observation vector then 


ote 
ee 
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This allows for truncation of the noise by removing the 
terms with small variance (i) as shown by Preisendorfer and 


Barnett (1977). 





GS 








WIT AWA 


WIT NWA 


sdola 


Sdold 


TOUd LA 


ANI LNOY YALAdWOO 


185 PReG pe ,eEWTxouddy 

SQUSTOTIJSON SUT], pe,Pouny, 
SUOJOIAUSBTY pazeounNI 

SLUSTOTJJIOD duvVTY, 

SUOJDIAUSBTY 

SONTeBAUSSTY 

XTALP_Y] SOUPTUPAOD OTALSWWAS PaZzTTeWUON 
1929S ePIeq SsnoTewoUuYy 


IWVN 


(HOS°OZT) V 
CHOS*OT). 


COL Wao) a 


CHOS 021) 
CO Ui) 
(OE sUcl) 
COC l= Gay 
(hOS*OZT) 


AVAYV 


O 


O 


a 


Al 


qd 


V 


-BeEP , 2eedjJ-SSstou, oy sonNpoud of 


iio tai eles UOLiOUnTs PEeUOcOUIIG™ PeOLATGhg soll Ust 1 atlooor vot clI> Ele p rot sao 


pesn sdom YeYR SoeuT Nod Jaynduoo/suotzedasdo jo souanbes ay} ST BuUTMOTTOF day], 


dq XITCNdddV 


66 





ody edeOBIL C1 





Figure Cl. 850 mb map types 1-6. MIsopleths are 
D-values (m). 


67 











LG 
ae 
Eh 
is 
ae 


SdsVO JO NOLLVIASC 
UTYWNN CYVCUNVIS 


Loe 
OT 
c 

2 

6 

5 
i 
ims 
Me 
8E 
ies 
62 
DOT 


sdSvVO dO NOTLLVIASd 
YdHWAN CYVAUNVIS 


“TPAdeUT YSPOsIOJ pue odAy dew Aq (°*TW°U) SUOTSTAAapP PTepueIS pue STOUTS qU QSg 


Geet 
SS a) 


Ie shore 
One 


Hh SO? 
Dore 
0°OOE 
el Ur 
8 la 
tee 
Suce 


uH-08 


Post 
Ca IGg 
O 39 
GC Ubr 
0°26 

8°26 

err le 
eae 
Sl (OIG 
8° 202 
Sook 
eo 
Ce ainell 


YH-hS 


0° é6h 
= Gshs 


8° 818 
oe WES, 


T° Sh 
8° 20S 
T°6ES 
E°OLS 
G°8Lh 
8° Shh 
6° onh 
NVdW 


€° 682 
4° En 
L£°n8S 
Calige 
Oe isk 
EG 
e608 
cb bene 
ST Shere 
€° 682 
a sinic 
OSC 
be JIC 
NVAW 


L¢ 


T 


0) 6O @) Wot 1. oe 


82 
ane 
Lh 
Ig 
62 


SdsVO JO NOLLVIAKd 
YdEWAN  CayVCNV.LS 


68h 
eal 
S 
OL 
Gil 
eal 
ile 
8L 
th 
LS 
68 
66 
OCL 


SdSvO JO NOTLIVIAAC 
CYVCNV.LS 


Ya eWAN 


C€° L87¢ 
DAs) 


Gee 
0° 962 


0° On? 
Gee 
i CL 
8°OcE 
L°8hzZ 
O°¢cL¢ 
8°8S2 


YH-82 


Oeil 
O'6dL 
OOS IE 
6°20L 
aes) 

Tes 

eats 

eelel 
ee OE 
G CUE 
L CGM 
ihe S100. 
Gmc 


lelmiOls 


tie Sill 
WN) aus 


Fame ihshe 
OmeOoG 


6°69h 
8° LES 
6°69h 
BGs 
Oy Guin 
le cin 
8° ETh 
NV 


legit 
1° BSZ 
O0°SSZ 
6°86 
SC onl 
Leste 
aoe 
6 GL¢ 
Se al 
c*TSl 
Oe L 
n° 9ST 
C°L9L 
NVAW 


de 


I 
il 
ee 
Ge 
Go 
Sie 
86 


ao as CO) ws Ae S&S 


SdSVO JO NOTIVIAG 
YdedWAN  CaVANVIS 


hOS 
oil 
L 
1 
HL 
ST 
ee 
GG 
oH 
9S 
[6 
LOL 
I 


SdsvVO JO NOLLVIAKG 
Ca VANV.LS 


UaeWNN 


LUG 
Ceca 


OVeGe 
syedilc: 

0 ine 
Pe ane 
le a All 
sh ele 
OS 8G 
OF Gee 
8° c8l 
8° S8T 


als) 018 


iS eke 
5°68 
0°69 
vas yf 
Oe 
i Od 
Bg 
O°6IT 
6° 26 
8° 08 
Boll 
Let 
9°66 


alls llke 


CROCE 
OP eaS 


Se oan 
mone 
Cech 
L° 6hE 
9° 00h 
ese 
LCE 
oi isre 
We te) ene 
a OUs 
NVAN 


6 Lou 
3° 98T 
IES 
6°SHTL 
2) SGA 
Le Gine 
8° LEL 
Gel 
ely Sheu 
Geol 
Os tone IE 
eel 
Creel 
NVA 


fal OOM Ceca olan CO at CO) 


ddAd 
dVW 


“TO TIAL 








7 








Fig. C2. cont. 700 mb map types 7-12. 
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Figure Dl. Examples of 500 mb maps. Isopleths are D-values 


CUM iC owe acc teteaeTOneby MAD type is andi— 
cated by the number above the map. 
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Pig. Dl. cont. Examples of map types 5-8. 
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